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Abstract 

We develop a cooling method controlled by a physical cooling radius that defines 
a scale below which fluctuations are smoothed out while leaving physics unchanged 
at all larger scales. We apply this method to study topological properties of lattice 
gauge theories, in particular the behaviour of instantons, dislocations and instanton- 
anti-instanton pairs. Monte Carlo results for the SU(2) topology are presented. We 
find that the method provides a means to prevent instanton-anti-instanton annihi- 
lation under cooling. While the instanton sizes are largely independent from the 
smoothing scale, the density and pair separations are determined by the particular 
choice made for this quantity. We discuss the questions this raises for the "physical- 
ity" of these concepts. 
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1 Introduction 



Beyond providing quantitative estimates for various physical parameters like masses and 
decay constants, Monte Carlo simulations of lattice gauge theories are increasingly involved 
in obtaining a physical picture of the structure of quantum field theories, such as the 
quantum vacuum and its topological properties. Such analyses are hampered by the 
ultraviolet fluctuations of the fields at scales comparable with the lattice spacing. A 
typical example of these problems is the scaling behaviour of the SU(iV = 2, 3) topological 
susceptibility which, for the commonly used lattice actions, turns out to be spoiled by the 
presence of topological fluctuations at the scale of the cut-off (dislocations). 

The need to improve the continuum approach of physical quantities has motivated 
the development of techniques to smooth out the UV, cut-off-dependent, structure of 
the configurations before measuring on them the quantities of interest, be these physical 
parameters like masses or features of the quantum vacuum. Apart from the improvement 
of lattice actions there are essentially two approaches that have been taken to realize this 
program: 

(i) preparing operators that are insensitive to UV fluctuations, and 

(ii) smoothing the fields themselves. 

For example, by smearing the spatial link variables [Q on fixed time slices one con- 
structs spatially extended operators with increased overlap onto the physical states, and 
thereby also achieves (i). This procedure preserves the transfer matrix and hence the phys- 
ical content. It has been successfully applied to potential and spectrum measurements. 

Smearing of the link variables can also be performed in an isotropic way, i.e. involving 
also the links in the time direction. When this is done on all links of the lattice before 
the measurements, it amounts to (ii). In this case locality and hence the transfer matrix 
interpretation are lost at the level of the original lattice. Isotropic smearing has been found 
to be useful to minimize UV renormalization effects of the charge operator in measurements 
of the topological susceptibility 0] and other properties of the gauge theory vacuum ||. 
However, it is not always clear if and to which limit the procedure converges [|[] . 

A special case of the second approach is cooling |J. It consists of an iterative, local 
minimization of the action, which proceeds sweep-wise and converges to configurations 
fulfilling the classical (Euclidean) equations of motion. In the case of continuum gauge 
theories, these are configurations characterized by the topological sector Q obeying S = 
\Q\Si, where S\ = j g 1 is the action of one instanton. Hence, up to scaling violations 
associated with the lattice discretization and assuming that topology preserves its meaning 
away from the semi-classical approximation, cooling should reveal the topological sector 
to which the original configuration belongs. 

Different smoothing techniques such as smearing H, cooling ||-[15] and inverse block- 
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ing [T9, O have also been applied in investigating whether gauge theory configurations 



can be described as instanton ensembles. Although results obtained for the topological 
susceptibility seem to be procedure- independent (i.e. they agree within 10%), the situa- 
tion is far less satisfactory for the determination of the instanton size distribution or the 
instanton density, which is mostly due to ambiguities in extracting the instanton content 



from the Monte Carlo configurations (for recent reviews, see [T 

Cooling is known to proceed as a diffusion process in the sense that the length scale 
up to which it smoothes the configurations grows with the number of iterations fl20| , |21 



There are two problematic features that are common to standard cooling procedures: 

1. In general the physical spectrum cannot be extracted from correlation distances 
shorter than the "cooling radius" to which smoothing has advanced. For example, 
the string tension extracted from a given pair of time slices rapidly diminishes as 
cooling proceeds (for discussions of this point, see below and [|TD|, 



Instanton-anti-instanton pairs, unlike superpositions of only instantons (I) or only 
anti-instantons (A), do not represent saddle points of the action. They belong to the 
topologically trivial sector and can have any action between 25*1 and 0. Consequently 
they are indistinguishable from trivial topological fluctuations and are removed by 
the cooling process. On the other hand, since in many instanton-based models they 
are conjectured to be relevant for important physical effects (for a review see P5| |) 



one would like to investigate if and how they are represented in the configurations. 

The question then is whether it is possible to apply a well defined amount of cooling, 
such that the configurations do not lose their physical properties at the length scales one is 
studying, but have become smooth enough to allow their identification. With the present 
cooling techniques, a well defined criterion to achieve this is lacking. The amount of 
cooling employed has to be fixed by subjectively judging when the configurations appear 
smooth enough to extract topological quantities^. In this paper, we introduce a systematic 
way of determining the cooling radius and relating it directly to a physical (continuum) 
length scale. This allows a smoothing of the fields up to a predefined scale, leaving all 
properties on larger scales untouched, and hence avoids the necessity of monitoring the 
cooling process, with all uncertainties involved in it. 

The procedure we shall introduce and discuss here is to let cooling saturate according 

to a local criterion. Following a suggestion by Niedermayer the local updating of links 

in the direction of minimization of the local action is only performed if the equations of 

motion are locally violated by more than a threshold 5 (see section 2). As a result, cooling 

1 We stress that this difficulty does not affect the topological susceptibility, measured on cooled con- 
figurations, since I-A annihilation does not change the topological sector Q. 
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will only proceed where the links are still far from the local minimum, and stop when the 
structure has everywhere approached the classical solution to a degree related to 5. As we 
will see, this allows to preserve the string tension above a distance given by the cooling 
radius and provides a natural lattice implementation of the valley method constraint to 
prevent instanton-anti-instanton annihilation (see e.g. [p5| ). 

In section 2 we introduce our cooling algorithm. The problem of relating the cooling 
radius to a continuum length scale will be discussed in section 3. In section 4 we test the 
effect of the algorithm on small instantons and I-A pairs, whereas section 5 contains results 
for the topological susceptibility and the properties of the instanton ensemble. Finally, 
section 6 is reserved for conclusions and outlook. 

The numerical results presented in this work have been obtained on the lattices in- 
dicated in Table |]. The 12 3 x 36 lattice has periodic boundary conditions in all four 



Lattice 


b.c. 


P 


a (fm) 


Confs. 


Sep. sweeps 


Therm, sweeps 


12 3 x 36 


Periodic 


2.4 


0.12 


800 


100 


20000 


24 4 


Twist in time 


2.6 


0.06 


350 


200 


20000 



Table 1: Simulation parameters 



directions, while the 24 4 lattice has twisted boundary conditions in the time direction, 
with twist k = (1, 1, 1). The large "time" extension and the twisted boundary conditions 
have been adopted in order to suppress finite-size effects peculiar to the \Q\ — 1 charge 
sector |37|, |T2"| . For the generation of configurations, we use the Wilson action. 



2 Improved cooling and restricted improved cooling 



In |p , 12| we introduced "improved cooling" (IC), based on a lattice action with improved 
scale invariance, as a well defined method to systematically eliminate UV noise and dislo- 
cations while preserving the topological charge associated with instantons above a certain 
size threshold po ~ 2.3a. Here we shall recall some results of IC that are relevant for our 
discussion. 

IC uses an action whose instanton solutions are scale-invariant for sizes beyond the 
dislocation threshold po of the order of the lattice spacing a. From a tree-level analysis pT 
one can determine a five-loop combination 

x+nu 

S m ,n = -^£ReTr(l- P~l ) 



4 



5 
i=l 

{mi, m) = (1, 1), (2, 2), (1, 2), (1, 3), (3, 3) for i = 1, . . . , 5 and: 

ci = (19-55c 5 )/9, pa = (1 - 64 c 5 )/9, 

c 3 = (-64 + 640 c 5 )/45, c 4 = 1/5 - 2 c 5 , (2) 

which has no tree-level 0(a 2 ) and 0(a 4 ) corrections at any 05. We have chosen C5 = 1/20, 
which appears to minimize the O (a 6 ) corrections for instanton configurations. 
The cooling algorithm is derived from the equations of motion 

U^W^xY - W^x)U^xy = , (3) 

where W is the sum of staples connected to the link Up(x) in Eq. ([!]). For SU(2), for 
instance, cooling amounts to the substitution (for clarity we drop the indices when possi- 
ble): 

U -> U' = V = W/\\W\\ , \\W\\ 2 = -Tr(WW*) . (4) 

Above the dislocation threshold, p ~ 2.3a, instantons are completely stable to any degree 
of cooling (at least if they are not affected by finite-size effects; see fl2| for a discussion 



of this point). The corresponding improved charge density using the same combination of 
loops leads to an integer charge already after a few cooling sweeps. 

We now supplement this cooling procedure by imposing the constraint that only those 
links be updated, which violate the equation of motion by more than some 5: 

U -> 17' = V if Aj(x) = ^ReTr(WW^ - (UW^) 2 ) > 5 2 , (5) 

where A M (x) is the dimensionful square of the lattice equations of motion. In the contin- 
uum limit (see f28|) 

Al(x) - -2 Tr((D u F Ufl (x)) 2 ) ; (6) 

hence we are locally requiring that the continuum equations of motion be satisfied to a 
degree specified by 5. This corresponds to a constrained minimization problem. For 5 = 
cooling proceeds unrestricted, i.e. all links are updated at every iteration, whereas for large 
values of 5 no changes are made at all and the configuration remains uncooled. Variation 
of 5 then admits a smooth interpolation between these extreme cases. 

In particular, for small 5 > the configuration is changed only locally, where it is 
far from a classical solution, but freezes where it has approached a solution to a certain 
approximation specified by 5 through Eq. ([Sp. Since the equations of motion contain 
gradients of the fields, 5 controls the energy of the fluctuations around classical solutions 
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and acts as a filter for short wavelengths. Hence cooling ends in a configuration slightly 
above the minimal action, the additional "action" being distributed smoothly over the 
lattice. We call this constrained variant of IC "restricted improved cooling" (RIC). 

In this paper we have used a slightly different version of (|5|), specific for SU(2): we 
replace fl5|) by the condition: 

U^U' = V if A M (x) 2 = \tt(1 - UV*) > 5 2 , (7) 



a 6 



with V as in (El). In the continuum limit 5 oc 5 



Lattice 


a 3 5 


5 

(fm" 3 ) 


Saturation 
at sweep No. 


a = 0.12 fm 


0.08 


46.30 


5 


n 


0.04 


23.15 


7 


n 


0.02 


11.57 


11 


n 


0.0085 


4.92 


22 


» 


0.005 


2.89 


40 


a = 0.06 fm 


0.04/8 


23.15 


20 


n 


0.02/8 


11.57 


35 


n 


0.085/8 


4.92 


« 85 


» 


0.005/8 


2.89 


« 180 



Table 2: The number of cooling sweeps after which saturation is reached for RIC. 



Since RIC does not update links that are sufficiently close to a solution, it changes 
fewer links after every iteration until no further changes occur. On the MC configurations 
we define saturation to be reached when less than 1% of the links are updated in a sweep 
and end the cooling thereQ. The number of cooling sweeps required to reach saturation is 
determined by 5. Some example values, for our SU(2) lattices, are given in Table ||. In 
Fig. |l] the fraction v(n) of links updated per cooling iteration is shown for the a = 0.06 fm 
lattice at two different values of 5. The figure elucidates the difference between IC stopped 
at a certain number of sweeps and RIC. For the former, of course, v is a step function 
whereas it is smeared out for RIC. 

2 We put this limit for practical reasons. We checked that it has little influence on the results: for 
very rough configurations (large 8) some narrow instantons may still shrink under cooling beyond this 
threshold; the effect is, however, at most 1-2% and disappears for smaller 8. 
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Figure 1: Fraction u(n) of links updated per RIC iteration as a function of n for the 
a = 0.06 fm lattice. 

3 Setting a scale for cooling 

In order to establish a physical scale up to which cooling is effective, we study the behaviour 
of the string tension calculated on the cooled configurations. Such analyses have been 
carried out for other variants of cooling in [pi], |22| , indicating that the string tension is 
lost over increasingly large distances while cooling proceeds. In RIC, however, the amount 
of cooling is related locally to the equations of motion, which have a continuum limit. 
Hence, we expect to be able to find a scaling criterion for tuning the cooling parameter S 
and thus define it directly in terms of a physical scale. 

A standard way to calculate the string tension is to extract it from the exponential 
decay of the correlation function of spatial Polyakov loops P^ L ' of length L according to 

@ 

^ (FW(x)pWt(O)) aM P (L)=a 2 a L L, aV 0O = aV i + ^. (8) 

X 

In order to see when the asymptotic value aMp(L) is reached, one may look for plateaux 
in the effective mass on a given time slice, 

aM(t) = -iln P ^ ( g ( + 1) . (9) 

For uncooled configurations, the signal is too noisy to extract the exponential fall-off. This 
problem is remedied by employing fuzzing techniques [0] to produce smeared operators 
with better overlaps onto the ground state. For the cooled configurations, this is usually 
not necessary as each cooling update includes the equivalent of one fuzzing step. 

In Fig. |2| we show, for the a = 0.12 fm lattice in Table |l|, the effective mass of the 
Polyakov loop for various numbers of cooling sweeps using IC or RIC (5 = 2.89 fm~ 3 ). 
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Figure 2: Behaviour of effective masses of the Polyakov loop under cooling. Left: IC. 
Right: RIC with 5 = 2.89 fm~ 3 . The dashed lines are to guide the eye. The horizontal 
line gives a string tension value (within the corresponding error band) as determined from 



Since in this exploratory study we have rather poor statistics, we do not try to obtain 
a precise value for the string tension but are only interested in comparing the plateaux 
obtained from the uncooled and cooled correlation functions. For illustration we indicate 
in the figure the bandwidth of the value obtained from other data sets for the string tension 
|30|j . The previously reported effect of a diminished string tension through cooling [ID , 22 ] 
is clearly visible, and in both cases there is a "cooling radius" at which the effective 
masses merge with their counterparts from uncooled configurations, and beyond which 
the string tension is still unaffected. However, in the case of IC (as well as other cooling 
variants), this cooling radius moves to larger distances with the square root of the number 
of iterations ||21[| . As a consequence, the physical properties of the system are changed up 
to distances that are so large that beyond them any measurement of physical observables 
with reasonable statistical errors is practically impossible. In contrast, RIC saturates in 
the sense that no more changes take place above a certain number of iterations. The 
physical properties of the system beyond the cooling radius reached at saturation appear 
to remain intact. 

Rather than by the number of iterations, the cooling radius for RIC is determined by 
the choice of the parameter 5, as displayed in Fig. [| A larger choice for 5 results in an 
"earlier" saturation, in the sense that links are left untouched when they represent a cruder 
approximation to the equations of motion, cf. Eq. @. As elucidated in the figure, this 
results in a smaller cooling radius. The clear approach to saturation offers the possibility 
to fix the desired cooling radius before a simulation and eliminates the need to monitor 
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Figure 3: The cooling radius of RIC as a function of 5. Left: 12 3 x 36, j3 = 2.4, 
a = 0.12 fin. Right: 24 4 , (3 = 2.6, a = 0.06 fin. 

the cooling process, which so far had to be stopped "by hand", based on a subjective 
judgement on the trade off between smoothness of configurations and conservation of 
physical properties. 

Furthermore, since the parameter 5 is expected to have a continuum limit, it should be 
possible to establish a systematic mapping between the values of S for a given lattice spac- 
ing and a physical length scale, up to which the smoothing process is active. This would 
permit a clean separation of physical scales into an ultraviolet sector, where smoothing is 
active, and an infrared sector, in which all physical properties are preserved. The naive 
expectation is that, close to the continuum limit, the dimensionful quantity 5 should scale, 
cf. eq. (!]). To check this behaviour we have repeated this study on the a = 0.06 fm lattice 
in Table p], which halves the lattice spacing of the previous one while preserving the same 
3-dimensional volume. The results for the effective mass on this lattice are presented in 
Fig. |^ (right). Despite the fact that the statistics is not very large, there is clear confir- 
mation of this expectation. We extract the cooling radius r(5) from Fig. |3] as the value of 
t at which the effective masses merge with the uncooled results and the string tension is 
recovered. A comparison of the data at fixed S for the two values of a shows (at the level 
of precision of these data) good scaling properties and therefore allows us to estimate the 
physical cooling radius to be 

r(S) « 0.8 <T 1/3 . (10) 

Since the string tension gives the fundamental physical scale of the theory, we con- 
jecture that the cooling radius defined in this way is valid for all physical observables, in 
the sense that all physical properties defined on scales larger than r(S) can be measured 



9 



on configurations Rl-cooled with 5. In the following we shall apply this conjecture to the 
observation of the topological structure. 

4 Instantons and RIC 

In order to provide a reliable extraction of the topological information from the MC con- 
figurations, there are two properties a smoothing algorithm should have: 

(1) It should allow a clean separation between dislocations (small instantons at the 
scale of the cut off) and physical instantons. The former should be eliminated under 
smoothing, the latter should remain unaltered. 

(2) Once we have fixed a physical distance below which smoothing is active, it should 
guarantee that the structure above it, among others instanton-anti- instanton pairs, re- 
mains unchanged. 

Improved cooling was designed to achieve (1) and we shall prove below that RIC 
preserves this property. With respect to the string tension, RIC also satisfies requirement 
(2). As we shall see this is also the case for I-A pairs. 

We present below the study of a set of prepared small instantons and I-A pairs that 
give an indication on the performance of RIC with respect to these two points. The reader 
interested only in the Monte Carlo results may skip this section and go directly to section 
5. 

4.1 Dislocations and RIC 

The improved action in Eqs. ([!]), @ was designed to minimize violations of scale in- 
variance: it is almost independent of the instanton size p when the latter is above the 
threshold po ~ 2.3a. Instantons with p > p are practically left unchanged under cooling, 
irrespective of whether we use IC or RIC. 

The situation is, however, quite different for small instantons with p < p$. Under IC 
these instantons shrink and disappear through the lattice "holes" . Such small instantons 
are expected to deviate considerably from solutions to the equations of motion. To check 
this conjecture we have evaluated this deviation for a set of instantons with size below Pq. 
We present on Fig. [| (left) the dependence on p/a of the lattice dimensionless quantity: 

Afat ak = yE^(A (11) 

with A M (x) from Eq. (|7|), evaluated on top of the instanton. As expected, for p < p 
the deviation from a solution increases in a rather steep way with decreasing instanton 
size. RIC limits A M (x) < 5, Vx; in the figure we indicate by horizontal dotted lines the 
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Figure 4: Left: Dimensionless deviation from a solution of the equations of motion: Af a ^ k 
- see Eq. {[HQ - as a function of the instanton size in lattice units. Right: Action of small 
instantons as a function of Af^ tt . The horizontal dashed line indicates the entropy bound: 
instantons with Swuson below this line should be considered as dislocations fi3l\] . On the left 
(right) figure the horizontal (vertical) dotted lines show the values of a 3 5 used in our RIC 
simulations. Open (filled) symbols correspond to the a = 0.12 (0.06) fm lattice. 



values of a?5 used in our RIC simulations. Instantons with Af a c t ak above each line will be 
destroyed under RIC with the corresponding 5. We may as well say that instantons with 
size below pQ IC (a 3 5) will shrink and disappear under RIC with 5, where Pq IC (o?5) is given 
by the intersections of the curves on Fig. |] (left) with the horizontal lines. 

We present on the figure data for both a = 0.12 and 0.06 fm. It is important to note 
that we are plotting the deviation from a solution evaluated in lattice (not physical) units. 
The rather weak dependence on a for p/a < 2 indicates that for instantons close to the 
lattice spacing the scaling law given by Eq. (^) does not hold, and the deviation from a 
solution is fixed solely in terms of the size in lattice units, i.e. p/a. 

To summarize, setting 5 allows to fix a threshold pQ IC (a 3 5) on the size of the instantons 
that will be preserved under RIC. 5 = corresponds to the IC threshold p = Po IC (^ = 
0) ~ 2.3a. For 5^0, the dislocation threshold (in lattice units) is smaller for larger lattice 
spacing. 

As pointed out before, the existence of this small instanton threshold is important to 
avoid the presence of dislocations after cooling. A semi-classical argument due to Pugh 
and Teper |3T| indicates that lattice instantons with action below the so-called entropy 



bound, Siattice(p) < 487r 2 /ll, dominate the path integral in the continuum limit and give 



rise to the lattice artefacts known as dislocations; for details of the argument, see [31 
Since we have generated the MC configurations with the Wilson action (Swuson) this is 
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the relevant quantity to characterize small instantons as dislocations. In Fig. |] (right) we 
show the improved action as well as Swnson f° r the instantons in Fig. |] (left). According 
to the above, instantons with 5wiison below the entropy bound (indicated on the figure 
by the dashed horizontal line) should be considered as dislocations. As can be seen from 
the figure the deviation from a solution of the equations of motion for such configurations 
is quite large. It is then obvious that they can be eliminated by putting a threshold on 
the maximally allowed deviation and this is precisely what RIC does. We indicate in the 
figure by dotted vertical lines the values of a 3 5 used in our RIC simulations. It seems that 
topological lattice artefacts are under control with RIC as long as we keep a 3 5 < 0.04. 



4.2 I-A pairs and RIC. 

One of the motivations for RIC was to introduce a physical scale above which the physical 
information concerning I-A pairs remained unaltered. Both on the lattice and in the 
continuum, these configurations are not solutions of the classical equations of motion, 
unless the instanton and anti-instanton are infinitely separated from each other. Since 
ordinary cooling is an unconstrained local minimization of the action, when cooled with 
it the I and A annihilate and decay into the trivial vacuum. The information concerning 
I-A pairs extracted from ordinary cooling is hence distorted in an uncontrolled way. RIC 
cooling, on the contrary, will leave unchanged those pairs for which the deviation from a 
solution remains everywhere smaller than 5. The expectation is that this provides a means 
to fix a physical scale above which the information concerning pairs is unchanged. 

To test this hypothesis we have analysed a prepared set of I-A pairs (similar studies 



for other variants of cooling have been performed in |fL5|j ). As in [15| we create a pair 
by putting together a Q — 1 instanton and, at a distance diA along the time direction, 
the anti-instanton obtained from it by time reversal. This configuration is I-cooled and 
its evolution is studied under cooling. We have monitored in particular the change of the 
deviation from a solution as evaluated at the instanton centre (x°) and at the mid-point 
between I and A (x m ): 

1 3 

A°' m = - M^), with x = x°, x m , (12) 
3 i=i 

where Aj(x) is given by Eq. (0) (the deviation measured on links along the separation 
axis is much smaller, therefore we did not include it in the estimate of A 0,m ). The results 
are presented in Fig. [5] (left). This figure should be interpreted as providing a series of 
I-A configurations labelled by the number of IC sweeps. As expected, with proceeding IC 
A°' m increases (the configuration gets distorted away from a solution) until it decays into 
a trivial one; after that we can no longer recognize I and A as peaks in the energy density, 
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Figure 5: Left: Evolution of the deviation from a solution of the equations of motion A°' m 
- see Eq. (jlW) - under IC as a function of the number of cooling sweeps. Right: A 0,m vs 



S(£ - see Eq. The results are obtained on the a = 0.12 fm lattice. The horizontal 

dotted lines show the values of 5 used in our simulations. Pairs for which A At (x) is below 
the line for all x are preserved under RIC with 5. 



the configuration approaches the trivial vacuum and A°' m decreases again. This picture 
changes under RIC. We have plotted in the figure horizontal dotted lines corresponding 
to the values of 5 used in our simulations. Whenever a configuration satisfies A M (x) < 5 
everywhere, in particular at x° and x m , it will be left untouched by RIC, i.e. the pair will 
be preserved and will not annihilate under RIC with 5. Pairs with A 0,m above the line 
will, however, be destroyed. 

Our objective is now to fix which pairs are preserved, using a physical length scale, 
since this may allow us to relate their effects to the physical quantities characterized by 
that scale. As already mentioned, since I-A pairs are continuously connected with the 
trivial vacuum, there is an ambiguity (in the continuum) in distinguishing a pair from a 
trivial fluctuation. The threshold separating the two is usually rather arbitrary and set 
"by hand" . One possibility is to do it through the value of the interaction between I and 
A, 

= 167T 2 - S IA . (13) 

It seems only justified to speak of an I-A pair instead of a fluctuation if Sj^ remains 
considerably smaller than the action of a "non- interacting" pair: I6n 2 . We can thus 
parametrize I-A pairs in terms of how small S}^ is, and use it as a threshold for which 
pairs will be considered as physical. We have monitored for the pairs in Fig. [| (left) 
and the results are presented in Fig. [| (right). Whenever we are able to still recognize the 
I and A there is a one to one correspondence between the distance to a solution (notice 
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A°' m is a dimensionful quantity) and Sj^. 

If we parametrize I-A pairs through how small Sj^ is, it is clear that setting 6 is 
equivalent to putting a cut on which pairs will be considered as physical and which will 
be thrown away according to the value set for the maximal Sj^-. It is in this sense that 
RIC implements the valley method constraint, since it amounts to a local constrained 
minimization of the action that prevents the interaction to act along the valley between I 
and A, if such an interaction is smaller than the chosen Sj^. 

5 Monte Carlo results 

5.1 Topological susceptibility from RIC 
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11.9(9) 


210.0(2) 


2.8(1) 


n 


0.02 


11.57 


199(3) 


13.2(8) 


109.7(2) 


2.93(9) 


n 


0.0085 


4.92 


198(2) 


13.0(6) 


51.7(1) 


2.91(6) 


n 


0.005 


2.89 


197(3) 


12.6(8) 


33.3(3) 


2.85(8) 


n 


IC nc=20 




199(2) 


13.1(6) 


25.2(2) 


2.91(8) 


» 


IC nc=40 




197(2) 


12.6(7) 


15.4(2) 


2.86(8) 


a = 0.06 fm 


no cooling 




297(4) 


22(1) 


50482(4) 


3.7(1) 


» 


0.04/8 


23.15 


203(7) 


4.7(6) 


80.4(2) 


1.7(1) 


!) 


0.02/8 


11.57 


203(7) 


4.8(5) 


40.3(2) 


1.7(1) 


J5 


0.0085/8 


4.92 


210(7) 


5.4(5) 


18.8(2) 


1.8(1) 


55 


IC nc=20 




205(7) 


4.9(7) 


30.8(2) 


1.7(1) 


J) 


IC nc=50 




205(7) 


5.0(7) 


12.9(2) 


1.7(1) 



Table 3: Topological susceptibility x, average charge squared {Q 2 ), action {S)/87i 2 and 
absolute value of the topological charge (\Q\), from IC and RIC. 

The results for the topological susceptibility are presented in Table |3| and Fig. ||]. For 
comparison we also include in the table and the figure the results for the topological suscep- 
tibility extracted directly from the MC configurations (without smoothing). As expected, 
the raw Monte Carlo data indicate strong deviations from scaling which, comparing to the 
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Figure 6: Topological susceptibility with IC and RIC. The results from RIC are plotted at 
the number of cooling sweeps corresponding to saturation - see Table [|. 



IC results, are particularly severe for the finer lattice 0. 

The expectation expressed in section 4.2 that changing a 3 5 allows to change the thresh- 
old on small instantons, and hence on dislocations, is corroborated by the MC data. The 
result on the a = 0.12 fm lattice at a 3 5 = 0.08 is clearly affected by dislocations, as seen in 
the tendency to approach the uncooled results. According to this observation and in view 
of the strong scaling violations of the uncooled data, it seems quite dangerous to extract 
the value of the topological susceptibility by extrapolating to zero smoothing. Generally, 
a minimal amount of smoothing may be indispensable to obtain a physically meaningful 
value for the susceptibility; see also fl32| . 

By contrast, in the "safe" region (a 3 5 < 0.04) we observe good agreement between 
the different results, including good scaling behaviour when comparing the a = 0.12 and 
0.06 fm lattices. This is an indication that most of the relevant physical instantons have 
been correctly taken into account and have sizes above the small instanton thresholds 
p* IC (a 3 5 < 0.04). 

The results agree well with the Witten-Veneziano formula and with results obtained 
previously by us using IC 0, as well as by other groups using different methods 
|2|, [14], jlS] . Differences with the results in || appear only for the largest (3 value and 
amount to about 10% of the value of x- 



3 This is to be expected for MC configurations generated with the Wilson action, since the simple 
semi-classical argument quoted above |}f| indicates that the contribution of dislocations to the Wilson 
action partition function diverges in the continuum limit. 
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5.2 Instanton ensemble from RIC 



A first impression of the effect of RIC is already given by the average value of the action 
presented in Table ||. Assuming that after some smoothing the action is mostly saturated 
by \Q\ = 1 instantons, the average action gives an estimate of the total number of instan- 
tons. It is clear that, under the above assumption, large S keeps more instantons, i.e. less 
pairs are annihilated. Moreover the scaling behaviour of 5 is confirmed within at most 
10% difference by comparing the average value of the action per unit volume between the 
a = 0.12 fm lattice at 5 = 23.15, 11.57, 4.92 fm" 3 : S/V = 16.28(2), 8.78(2), 4.57(1) fm" 4 ; 
and the corresponding values on the a = 0.06 fm lattice: S/V = 18.49(4), 9.27(4), 4.32(4) 
fm -4 . 

To associate the action with an instanton ensemble one has to extract in some way 
the instanton content of the configurations. This is a difficult task due to several reasons: 

(i) as pointed out before, there is an ambiguity in distinguishing between an I-A pair 
and a trivial fluctuation; 

(ii) this ambiguity is enhanced if the ensemble is dense, as seems to be the case, since 
then the objects in the ensemble are strongly overlapping. 

This last point sheds quite some doubts on the reliability and usefulness of the in- 
stanton parametrization, unless one can justify the introduction of some physical criterion 
to select the relevant pairs for the process under consideration. In particular, the results 
obtained for instance for the size distribution seem to be quite dependent on the particular 
method used to smooth and analyse the configurations (see 



Our particular way of extracting the instanton information out of the configurations 
relies on two assumptions: (1) instantons should appear as local self-dual peaks in the 
action and charge density; (2) only pairs with Sf^ considerably smaller than 167r 2 should be 
considered as such. Hence we approximate the action and charge density by a superposition 
of self-dual non-interacting instantons and anti-instantons parametrized through the 1- 
instanton BPST ansatz (further details on the procedure are provided in the Appendix). 
We can measure the departure of the real action and charge density from the above non- 
interacting ansatz through the quantities 



Jd 4 x\q(x) -q ftt (x)\ 2 
\ f#x\q{z)\* { ] 



and e s , defined as e q with the replacement q — > s. In previous work we have adopted 
the prescription of performing a minimal amount of cooling such that the deviations from 
the ansatz described above were at most 0.3. Under this restriction, and using IC, we 
observed: 

(i) for the minimal level of smoothing, a rather dense ensemble with a number of 
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instantons by about a factor of 2 larger than the dilute gas prediction (N) = (Q 2 ), with 
iV the number of instantons plus anti-instantons; 

(ii) a size distribution peaked at about 0.4 fm and rather stable under further cooling (a 
much smaller average size of ~ 0.2 fm has been reported in [[3], |I6|, extracted from smeared 



configurations and under extrapolation to no-smearing; see 5.2.2 for further discussion); 

(iii) a rather homogeneous distribution of I and A, which seem to be randomly dis- 
tributed over the lattice. 

For the largest values of 5 in our RIC analysis we will present results corresponding to 
values of e s ^ q up to 0.6, where only a small part of the action and charge density is correctly 
described by our instanton ansatz. These results should be taken with care, since they 
correspond to rather dense ensembles and, as mentioned above, in such case there is a 
significant ambiguity in the extraction of the instanton content. The aim in the present 
work is to study the results obtained by our instanton finder as a function of 6 in RIC. 



5.2.1 Results 

Now we turn to the analysis of the Monte Carlo configurations. Results are presented in 
Table |. 



a 


5 (fm- 3 ) 


r(5) 


(N/V) 


(On) 


(Oia) 


(d n ) 


(diA) 




(p) 


0.12 


23.15 


0.24(5) 


14.86(7) 


0.95 


1.28 


0.39 


0.28 


0.75(0.55) 


0.38 


j) 


11.57 


0.35(5) 


8.72(2) 


0.86 


1.08 


0.46 


0.36 


0.60(0.45) 


0.40 


)) 


4.92 


0.50(5) 


4.21(2) 


0.79 


0.88 


0.54 


0.49 


0.45(0.35) 


0.42 


)) 


2.89 


0.60(5) 


2.70(3) 


0.76 


0.79 


0.59 


0.58 


0.35(0.30) 


0.44 


)) 


IC(20) 




2.06(2) 


0.70 


0.72 


0.61 


0.61 


0.30(0.30) 


0.42 


0.06 


23.15 


0.27(3) 


18.7(1) 


0.93 


1.20 


0.40 


0.30 


0.60(0.45) 


0.35 


)) 


11.57 


0.35(3) 


9.6(1) 


0.84 


1.02 


0.48 


0.39 


0.45(0.40) 


0.38 


n 


4.92 


0.55(3) 


4.35(6) 


0.77 


0.85 


0.57 


0.53 


0.25(0.25) 


0.40 


n 


IC(50) 




2.96(5) 


0.68 


0.71 


0.59 


0.57 


0.25(0.25) 


0.37 



Table 4: Cooling results for IC at 20 and 50 cooling sweeps and for RIC at various 5. 
The lattice spacing (a), the cooling radius (r(5)), distances and sizes are given in fm. The 
density (N/V) is in fm~ A . The errors are statistical, when they are not explicitly quoted 
they amount to 1 or less in the last indicated digit. 



We are mostly interested in studying the performance of RIC with respect to I-A pairs. 
First we observe that the density of instantons (N/V in Table |j) increases with 5. This 
corroborates the finding that less I-A annihilation takes place for large 5. Moreover the 
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x IC a = 0.06 fm 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 



p(fm) p(fm) 

Figure 7: Size distributions on the a = 0.12 fm lattice (squares, triangles) and a = 0.06 
fm lattice (crosses). The vertical dotted (dashed) lines indicate the location of the IC 
dislocation threshold po — 2.3a on the a = 0.12 (a = 0.06) fm lattice. 



density scales correctly (up to at most 10% (20% for the largest 5) difference) by comparing 
the results for the two lattice spacings at fixed S. 

We define pairs in this ensemble by assigning to each instanton the anti-instanton that 
has maximal overlap Oia with it; Oia is defined as: 

n Pl + PA , 

0lA = ^ovT' (15) 

where pi t A are the sizes and d/A the distance between I and A. The results for the average 
overlap and separation are given in Table |j. As expected, large 5 preserves stronger 
overlapping pairs. Notice the good agreement between the pair separation djA and the 
cooling radius r(5), Eq. (fTOj) , a clear indication that smoothing has been active in removing 
pairs with a separation smaller than the cooling radius. As we will further discuss below, 
this also indicates that there does not seem to be an intrinsic scale for I-A separation apart 
from the one given by r(5) itself. 
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For comparison we also present in Table |] the overlap and separation between maxi- 
mally overlapping objects of the same charge. Note that for large 5, i.e. denser ensemble, 
objects with opposite charges are closer on average than objects with equal charges. How- 
ever, in our earlier IC analysis [fT2 |, which roughly corresponds to RIC at small 5, we 
observed a homogeneous distribution of I and A over the lattice. The difference can be 
understood by noticing that pairs with strong overlap will be removed faster during cool- 
ing, and have probably been eliminated at the stage of smoothing provided by IC. Since 
the action of an I-A pair decreases as the overlap increases and is smaller than that of an 
I-I or an A-A pair, it is indeed energetically favourable that I-A pairs overlap more. 

In Fig. [7| we present the instanton size distribution as a function of 5, Table [| provides 
the average instanton size (p). Since IC and RIC converge to scale-invariant instanton 
configurations, the instanton sizes, in contrast to other physical quantities, should not 
change when the cooling radius increases beyond them (up to modifications that would be 
induced by I-A annihilation and possible dressing of the instantons by fluctuations). This 
is certainly the case for the results obtained with IC above the minimal smoothing level 
(see ||12|| ) and also for RIC with small 5 (i.e. density of instantons below 4 fm~ 4 ). Note, 
however, that now we go to much denser and rougher ensembles. At these densities the 
stability of the size distribution is not as apparent, both the average size and the shape 
of the distribution weakly depend on the smoothing level. The increase in the close pair 
population seems to lead to a certain shift in the average size, both by spoiling the size 
estimate and by re-weighting the size distributionF]. It is difficult to evaluate which of these 
two effects dominates, cf. the debate about this problem in the literature (for reviews, see 



18| , |19[). According to || smoothing the configurations induces a modification of the size 



and the correct one should be obtained by linear extrapolation to zero smoothing levels 
(for SU(3) this seems to give rough agreement between the results from different smoothing 
techniques |19|)- Our attitude is not to trust the size determination corresponding to large 
densities (small cooling radius), since, as we have discussed, it is difficult in such cases to 
obtain a reliable description of the action and charge density in terms of instantons (the 
deviations amount to 0.6 (0.5) of the total charge (action)). 

This problem cannot be easily solved and is indeed farther reaching, since it is related 
to the general question of whether there exists an effective I-A density relevant for a 
given physical effect. The numerical agreement between the cooling radius and the I- 
A separation - see Table § - suggests that the latter is determined by the former, i.e. 
there appears to be no intrinsic scale for I-A separation independent of the one set by the 
smoothing scale itself. 



4 Note that the size distributions of Fig. |t] are normalized, while describing instanton populations 
varying by a factor up to 10. 



19 



0.6 



0.2 



(p) 


- <d IA ) /: 


lii 






0.2 0.4 0.6 
r(5) fm 



0.2 0.4 0. 
r(6) fm 



- 20 



15 



10 



0.2 0.4 0.6 
r(<5) fm 



Figure 8: Dependence of the average size (triangles), the I-A average distance (dia- 
monds) and the instanton density (squares) on the cooling radius r. The dashed curve fits 
the density data to (N/V) = (A r /V r ) asymp totic + C/r 4 . Open (filled) symbols correspond to 
the a = 0.12 (0.06) fm lattice. 



In order to elaborate on this, we consider again the instanton properties as a function 
of the cooling radius. It is apparent from Fig. [7| and Fig. |8| that the dependence of the 
average size upon r(S) is rather weak, even at small r(S) (large 5). This seems to indicate 
the existence of a non-vanishing, dominant instanton size even at cooling radius zero. 
However, Fig. [8] also shows the different behaviour of the density (N/V) and of the pair 
separation (djA). A simple extrapolation of the data to zero cooling radius is compatible 
with a diverging density and a vanishing I-A distance. This is in accord with the fact 
that there is no principal distinction between a trivial fluctuation and an instanton-anti- 
instanton pair, apart from the somewhat arbitrary criterion given in terms of S 1 ^. This 
supports the observation made above that there appears to be no average I-A separation 
or instanton density that could be defined independently of the smoothing scale (at least 
up to the scales we have probed). Thus, speaking of an instanton ensemble seems to be 
possible only in the context of smoothing and by providing a specific smoothing scale. If 
this is indeed the case, a criterion for selecting "physically relevant" pairs may only come 
from considering the relevant scale of the phenomena to which we can expect these pairs 
to contribute. 



6 Conclusions and Outlook 

By fixing the amount of cooling with the help of a cooling radius (r(S)) above which phys- 
ical observables such as the string tension remain untouched, restricted improved cooling 
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provides a smoothing scale for treating MC configurations. Because of its connection with 
the equations of motion this scale appears as a continuum length scale, which can be 
systematically determined from the measurement of one physical quantity. Starting from 
the rough, UV-dominated field configurations (cooling radius r = 0), RIC with increasing 
r acts as a "low-pass filter" by smoothing out wavelengths smaller than r and producing 
configurations on which correlation functions at distances beyond r retain the physics of 
the larger wavelengths. The conjecture that the cooling radius r(S) provides a universal 
smoothing scale seems well supported by our analysis. 

The observation of topological structure itself poses a difficult problem. The topo- 
logical susceptibility appears to be quantitatively under control, but the more detailed 
topological structure raises many questions. Firstly, for very small cooling radius the ob- 
served structure is very dense and only poorly described in semi-classical terms. Secondly, 
we find that cooling uncovers topological structure on scales close to the cooling radius. 
While the instanton size distributions are rather stable under variations of the smoothing 
scale, the instanton density and the mean separation of instanton-anti-instanton pairs ap- 
pear to be determined by it. Decreasing the smoothing scale, this structure thus becomes 
denser until it continuously disappears into the UV fluctuations. Hence, for these quanti- 
ties, an extrapolation to no smoothing does not lead to meaningful results. On the other 
hand, it appears to be quite arbitrary to choose from the observed structure a "physical" 
sub-ensemble by fixing a finite smoothing scale. It may well be that the characteristic 
length scale of the phenomena under scrutiny has to be taken into account when estab- 
lishing the criteria for finding an "effective" relevant structure that contributes to these 
phenomena. Here RIC may provide a useful tool, since it introduces a smoothing radius 
as a continuum length scale, which can be systematically varied in the analysis. 

Further studies should provide improved statistics for the precise determination of the 
cooling radius and the other properties of RIC, but also a better theoretical understanding 
of these questions. 
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Appendix 

We provide a few details about the way we perform the extraction of the instanton content 
from the configurations (a more detailed description of the method will be presented 
elsewhere) . 

As already mentioned in section 5, our instanton identifier relies on the following 
assumptions: 

1. Instantons are self-dual objects with energy density localized in space and time. 
Self-duality imposes, in particular, that the localization of the peaks should be the 
same in the electric and magnetic parts of the action density and also in the charge 
density (for details on how we imposed the cut from self-duality, see [0). 

2. We approximate the action and charge density by a superposition of self-dual non- 
interacting instantons and anti-instantons parametrized through the analytic expres- 
sion for the 1-instanton action and charge density: 



Si(x) 



48 



i 



i + E 



(16) 



S &t(x) = S i( X ) 5 <7fit( X ) = r-J 2 n i S i( X )> ( 17 ) 

i=l,N ° n i=l,N 

with N = Nj + Na the number of instantons plus anti-instantons, p l and respec- 
tively the size and location of the ith (anti)instanton and n, = 1 for an instanton, 
rii = — 1 for an anti-instantonF]. 

The departure of the real action and charge from the above formula is measured 
through the quantities e q ^ s given in Eq. fll4|). For the analysis in this paper we have 
supplemented the one performed in Jl2]| with a further restriction: every time an instanton 



candidate is located, it is counted only if, by adding it to sg t and q^ tl s s , q simultaneously 
decrease. Once the instanton has been accepted, its contribution (sj(x)) is subtracted 
from the total density before proceeding to analyse the next instanton candidate. Note 
that this procedure assumes that all self-dual objects in the ensemble can be described in 
terms of a superposition of instantons following the BPST ansatz (other possible self-dual 
objects, which are not well parametrized by this ansatz, would be discarded by our e S)g 

5 To account for periodicity effects we add to the energy density of each instanton the eight closest 
replicas, obtained by displacing (M) by a torus period in one of the four directions. 
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cuts). A small value of e s>q is, hence, a signal that most of the structure present in the 
ensemble is well described in these terms. 

As an example of the kind of results we expect, let us consider the I-A pairs in Fig. ||, 
where e s (e q ) varies from about 0.1 (0.05), for the pair with smallest in the figure, to 
0.32 (0.18), for the one with largest Sj^. Beyond the latter, our instanton finder fails to 
identify I and A any longer. In the same way, for the small instantons in Fig. |4], e s (e q ) 
goes from 0.05 (0.15) for the largest one in the figure to 0.06 (0.86) for the smallest one. 
Again beyond that point our finder no longer sees an instanton at this location. 
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